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Abstract The problem of formation of the aerosol of space origin in Earth’s atmosphere is examined. 
Meteoroids of the mass range of 10' -10' g are considered as a source of its origin. The lower bound of 
the mass range is chosen according to the data presented in literature, the upper bound is determined in 
accordance with the theory of Whipple’s micrometeorites. Basing on the classical equations of 
deceleration and heating for small meteor bodies we have determined the maximal temperatures of the 
particles, and altitudes at which they reach critically low velocities, which can be called as “velocities of 
stopping”. As a condition for the transfonnation of a space particle into an aerosol one we have used the 
condition of non-reaching melting temperature of the meteoroid. The simplified equation of deceleration 
without earth gravity and barometric formula for the atmosphere density are used. In the equation of 
heat balance the energy loss for heating is neglected. The analytical solution of the simplified equations 
is used for the analysis. 

As an input parameter we have used the cumulative distribution of space matter influx onto earth 
on masses in large mass range. Basing on this distribution we have plotted three-dimensional probability 
density distribution of influx of particles as a function of parameters, which detennine the heating and 
stop altitude of a meteoroid: initial mass mo, velocity of entry into the atmosphere Do and radiant zenith 
angles zro. The obtained three-dimensional distribution had been presented first as a product of three 
independent distributions on the mentioned parameters, then it was transformed using the equation of 
deceleration into the distribution on the following parameters: mo, vo and “altitude of stopping” I Is. The 
final 2-dimensional distribution on parameters vo and //$ of the aerosols of space origin in the 
atmosphere was obtained by means of integration of the previous distribution over Do- 
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1 Introduction 

There are aerosols of both ground and space origin in Earth’s atmosphere. Aerosols of the ground origin 
are presented basically in the lower atmosphere: in the troposphere. The most powerful aerosol layer of 
the ground-based origin, known also as Junge Layer, is placed at altitudes of 10-25 km. It originated 
from the condensation of some components of the atmosphere appearing from the photo-chemical 
transfonnations of some products of volcano eruptions, for example sulphuric acid vapors. The second 
confidently established aerosol layer in the atmosphere is placed at altitudes of 80-85 km, corresponding 
to the minimal atmospheric temperature, in the mesopause. The origin of this aerosol layer in not finally 
established. Most of scientists, and the authors as well, hold an opinion that all the particles there to be 
of space origin. Under some special conditions the condensation of water vapors on these particles 
becomes possible, and we can see, probably, the high-latitudinal silvery clouds. 
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According to meteor physics investigations (Whipple 1950; Whipple 1951; Levin 1956; Opik 
1956; Lebedinets 1980; Lebedinets 1981; Voloshchuk et al. 1989) most of low-mass particles coming 
into the atmosphere with initial velocities -11.2-72.5 km /s lose their energy at altitudes of 140-80 km. 
Small fragments detaching from already heated bigger particles in the atmosphere cannot be decelerated 
without almost entire loss of their masses due to evaporation. The deeper penetration of a particle to the 
atmosphere the lower the probability to save its macro size. This task of motion, deceleration and 
destruction of a separated particle in “abnormal environment” according to the terminology of Opik 
(1956) we were considering in Voloshchuk et al. (1989). Such a conclusion is also given from the 
experimental investigations of chemical analysis of particles, caught in the atmosphere with the help of 
high airplanes and balloons. Such particles are similar to coaly chondritics (Nady 1975), having a big 
amount of helium in their surfaces, which penetrated there from the solar wind. Therefore, these are the 
primary interplanetary particles, which have come though the atmosphere without intensive heating and 
are not the products of fragmentation of larger bodies (Brownlee and Hodge 1973). 

The amount and distribution of the aerosol of space origin in the atmosphere is connected by 
some authors with planetary global warming. In this work we will try to examine the problem of 
formation of the aerosol of space origin in Earth’s atmosphere basing on the initial meteoroid 
distributions on the Earth’s heliocentric orbit and the equations of classic meteor physics. 


2 Meteor Physics Equations to be Used 

In this chapter we consider the basic equations of meteor physics to be used in the work, namely the 
equation of heating of the meteoroid, and the equation of its deceleration. In addition, the simplification 
of the equations in order to realize the final investigation analytically is substantiated. 

2.1 Complete Equations of Meteoroid Deceleration and Heating 

The base assumption for the transformation of a small meteoroid into an aerosol particle, not into a 
meteor, consists in non-reaching by the meteoroid its melting temperature. Therefore, we have to 
determine the mass interval, and other parameters of meteoroids, which coming into the Earth’s 
atmosphere, do not reach the melting temperature because of their deceleration and heat radiation. 

2.1.1 Heating Balance Equation 

The theory of heating of low-mass meteoroids with their deceleration, which plays an important role in 
this case, were developed by Whipple (1950), Whipple (1951) and later by Fecenkov (1955). They have 
obtained the name of Whipple’s micro-meteorites. It is known (Levin 1956) that the particles having the 
size less than xo warm up to the same temperature (xo is the wanning up depth at which the temperature 
of the body is less to e times relatively the surface). According to Opik (1937) and Levin (1956) such 
particles have radius ro < 10' cm. The change of temperature of such a particle with taking into account 
the energy loss for heating and radiation can be written as: 

S U0 Edt = m 0 cdT + fla ( T 4 - 7 0 4 )S FO dt, ( 1 ) 

where Smo = const and Sfo = const are the middle section and entire surface area of the particle 
accordingly, mo = const is its initial mass, c is the specific heat capacity and cr is Stefan’s constant, T and 
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To are the current temperature for time t and initial temperature of the particle in the field of solar 
radiation at the distance of 1 a.u., /? < 1 is a coefficient of thermal radiation of the meteoroid 
characterizing the digression from black body radiation, E = A p A v 12 is the energy incoming to unity of 
the meteoroid surface due to its collision with atmosphere molecules, A is the dimensionless coefficient 
of heat conductivity, p A is the atmosphere density. 

2.1.2 Deceleration Equation 

If the space particle is not warmed up to the melting temperature it becomes the aerosol particle. So, the 
next question we should answer: at which altitude will it stop? In order to solve this problem we 
consider the equation of deceleration, which can be written in the most common vector view as 

= - c r S moPaW + m (} g , (2) 

at 

or separated into constituent parts: 

m 0 -37 = -c R S MQ p A v 2 + m Q g cos (3) 

at 


v 


d Z R 

dt 


“S' sin z R , 


( 4 ) 


where v is the meteoroid velocity, zr is the zenith angle of its radiant, cr is the resistance coefficient, g is 
the free fall deceleration constant. 


2.2 Accepted Assumptions 

In our calculations we use some assumptions and simplifications. First, we suppose the particles of 
space origin producing the aerosol are of meteor mass range. Lower bound of meteoroid initial mass is 

jo 

10“ g according to meteoroid mass distributions presented in literature (Ceplecha et al. 1992), higher 
bound corresponds to the r 0 , and is approximately equal to 10' g (Opik 1937). The second assumption is 
that we consider just wanned up and evaporated particles and neglect the mass loss due to blowing 
meteoroid molecules away in its “cold” state. The next, the most doubtful assumption consists in the fact 
we use the barometric formula for the atmosphere density: 

P A (H) = p A (0) exp(- ^—) . (5) 

ti 

^ • * 

Here p A ( 0), 77 are the atmosphere density at the sea level and altitude of the homogeneous atmosphere 

accordingly. For precise calculations one should use the numerical solution of the equations (1) and (2) 

and take the real atmosphere density distribution from modem models of atmosphere, especially for 

altitudes over approximately 120 km. We use the formula (5) here just for the purpose of obtaining the 

analytical solution of (3) and (4) in order to understand the physics of the aerosol layer formation. Then, 
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we consider the sporadic meteoroids as the main source of aerosol particles, i.e. the particles which are 
supposed to be of the stone composition. Finally, we calculate the mean aerosol influx during a year. 

2.3 Simplification of the Equations 

According to Opik (1937), the small meteoroid spends almost all its energy for the thermal radiation if 
its radius r < 10' cm (corresponds to mo ~ 10" g for spherical particles), so we can neglect the first term 
in the equation (1): 


T 4 



Smo*Pa»- 

2/3ctS f0 


( 6 ) 


Since we deal with low-mass particles we can suppose they are decelerated relatively fast, so we 
can neglect the gravity tenn in the equations (2). The equations (3) and (4) in this case transform into the 
equation 


do _ 2 

m o ^ CrSmPaO ■ 


(7) 


and z R = z R0 = const . 

Also we use the relation between the time t and altitude H of the particle 

dH = -o cos z R dt (8) 


and express the middle section and surface area of the particle through the shape parameter A: A = 
Sm/V , where V is the meteoroid volume. Supposing the particle is spherical .SVo = 45 mo = 4A(m 0 /pM) , 

9 /9 

the shape parameter for spherical particles to be A = 7r(3/47r)' . 

2.4 Variation Parameters, Constants and Final Equations 


Using (6), (7), (8), the shape parameter and barometric formula (5) we obtain 

( 1 ATT* \ 


7; 4 =^^exp 


8/fo- 


3 c r AH 


V 


1/3 2/3 

m 0 p M cos z R0 


Pa 


(9) 


7 


f 


o = u 0 exp 


c r AH 


1/3 2/3 


v m o Pm cos z r 


~Pa 


( 10 ) 


Reaching by the particle of maximal temperature along its trajectory can be derived from dT/dp A = 0 , 
and so from (9): 


Pat 


max 


1/3 2/3 

m o Pm c °SZ 
3 c r AH* 
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Putting it back into (9) we obtain 


T 4 _ T 4 — Pm c °SZ R 0 v 0 

max 0 “ 24 j3ac R AH* exp(l) ' 


(11) 


Thus, the condition of transfonnation of space particle into the aerosol can be expressed now as 


T < T 

max melt ’ 


( 12 ) 


where T max has to be expressed from (11), T me / t is the melting temperature of the particle. 

Looking at the equations (10) and (1 1) we can note that there are three parameters of a meteoroid 
(under the assumptions made above) having an influence onto its belonging to the class of aerosols or 
meteors, and to the altitude of stopping in the case of the aerosols. These are initial mass of the particle 
mo, velocity vo, and zenith radiant angle zrq. The ranges of their variations are: mo = 10" - 10' g 

according to Ceplecha et al. (1992) for the lower limit and Opik (1937) for the higher limit (see above), 
vo = 1 1.2 - 72.5 km/s, i.e. the particles belonging to the solar system are considered, z R0 = 0° - 90°, all 
possible entrance angles are taken into account. 

Expressing (12) through the variation parameters we obtain the final inequality of separation of 
the meteoroids onto aerosols and meteors 


m^ol cosz R0 < C T , 


(13) 


where C r = 24Pctc r AH* exp(l)(C, -T Q 4 )/Ap™ . 

If the condition (13) is realized we can find the altitude of stopping Hs of the aerosol particle 
from (10), supposing the velocity of stopping w s is a small enough value. Here we continue to use the 
equation (10) except the Stokes fonnula for low velocities, where the deceleration is proportional to the 
first power of the velocity, so H$ can be found from the expression 


C v p A { O)exp(-^f) 


= v s exp 


H 


m T COSZ R0 


(14) 


where C v = c R AH*pJ n . 

During the calculations the following values of constants are taken (Levin, 1956): A = 1, o = 
5. 67032x10° erg-cm'^K'^s' 1 , p = 1, Cr = 1, H* = 7x 10 s cm, p^(0) = 1.6><10' 3 g/cm 3 , p M = 3 g/cm 3 , To = 
276 K, Tmeit = 1600 K, vs = 0.5 km/s. For an iron particle A = 0.75; pm = 7.6 g-cm°; c R = 1.25; T me u = 
1800 K. 


3 The Statistical Approach to the Process of Space Origin Aerosol Formation in the Atmosphere 

Here we propose the statistical model for the description of atmospheric aerosol formation from 
meteoroids. We will construct the 3-dimensional distribution on variation parameters having an 
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influence onto the probability of the aerosol formation and onto the altitude of the aerosol layer. Let’s 
represent the distribution as a multiplication of three independent single-parameter probability density 
distributions: 


Pmuz K > U 0’ Z Ro) = Pm K )Pu («0 )Pz ( Z RO ) • 


(15) 


It is obvious that this distribution should be normalized to unity over all three parameters, i.e. there must 
be 


J J \P m uz( m 0 ’°0’ Z R0) dm 0 dV 0 dZ R0 = 1 > 

where the integration is carried out inside all possible ranges of parameter values. 

3 . 1 Primary Distribution of Meteoroids 

Let us find all three 1 -dimensional primary probability density distributions, and start from the 
distribution on mass. 

3.1.1 Probability Density Distribution on Initial Mass 

There can be found in the literature distributions of space matter onto Earth as cumulative distributions 
of number of particles on their masses, for example Ceplecha (1992), Kruchynenko (2002), 
Kruchynenko (2004). We will use the linear dependence (Kruchynenko 2002, Kruchynenko 2004) for 
the further calculations: 

log 10 N(m 0 >m 0 ) = C 0 -k log 10 m 0 , (16) 


where N(md > in 0 ) is a number of particles with masses not less than m t) coming into all Earth 
atmosphere during a year, Co = 7.86, k = 0.892. 

The probability density distribution on mass p m (m o) according to cumulative distribution (16) can 
be described by Pareto distribution: 


P m K < m J = 0 

/ x knC, 

p m (fn 0 >m o/) = — jnr 

m o 


(17) 


where m 0 i is chosen freely. The probability density function is normalized to unity in the value range 0 - 
+oo. There are the following obvious consequences: 


0 fTl 

F( m o) = f Pm ( OT 0 ) dm 0 = 1 - f Pm K ) dm 0 = 1 f 


(18) 
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dN(m 0 ) 

N,(m 0 i<m 0 <+cc) 


dF (m 0 ) = p m (m 0 )dm 0 , 


(19) 


where F(/«o) is the cumulative probability, Nfano > /??o/) is a sample of all particles in the chosen range, 
which can be found from (16) as N,(m 0l ) = \0 C ° / m k ol . We suppose m 0 i = 10' 18 g, so N/ = 8.24><10 23 . 

3.1.2 Probability Density Distribution on Initial Velocity 

The probability density distribution on velocity p v (v o) will be chosen according to radar meteor 
observations, for example (Voloshchuk et al. 1989): 

Pu K ) = PG(Pi , ct o1 ) + (1 - P)G(U 2 , cr u2 ) , (20) 


where 


G(o) 


1 



exp(- 


(o-o ) 2 


are Gaussians with the following parameters: v 1 — 32.32 km/s, o v \ = 6.51 km/s, uj = 54.26 km/s, o V 2 = 
5.15 km/s. The value P is changing during a year. For the mean value we choose P ~ 0.33 (Voloshchuk 
et al. 1989). It is obvious that the probability density function is normalized to unity in the range 0 - +oo. 

3.1.3 Probability Density Distribution on Initial Radiant Zenith Angle 

The probability density distribution on radiant zenith angle p z (z R0 ) will be derived from the following 
thoughts: let suppose that the number of particles dN(r, r + dr) entering into earth atmosphere from some 
direction in the range dr in some spatial angle dQ (see Figure 1) per time unity can be expressed as dN(r, 
r + dr) ~ 2ni)7irdrdQdt, where no is a spatial concentration of meteoroids. Since r = /?© sin zr, we have 
dN(zR, z R + dz R ) ~ 2no7rR® 2 sinz R cosz R dz R dQdt. So we have to use the sine-cosine distribution sin z /m cos 
zro. 

After normalization to unity we obtain the final distribution on zenith radiant angle 
Pz(. z ro) ~ ^ sin z R0 cos z R0 . (21) 


Strictly saying, this distribution will be distorted by the Earth gravity, but we use it due to its simplicity. 



Figure 1 . To the derivation of the probability density function on radiant zenith angle. R© is Earth’s radius. 
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3.2 Separation of the Distribution into Aerosols and Meteors 


The primary distribution of meteoroids can be conceived as a geometrical 3-dimensional cube, where 
the three considered parameters mo, vo and z R o determine the dimensions along its three ribs-axes, limited 
by the permissible parameter values. The “intensity” in each point inside such a cube is expressed by the 
value of p, m): {m 0 , vo, z R 0 ). The real number of particles dN in the range of dm 0 , dv o, and dz R 0 can be found 
from (19). If we make a few sections perpendicularly to the cube rib describing the mass mo we will 
obtain the 2-dimensional pictures in coordinates vo <-> z R o for the fixed mass values, where the relative 
value of p mDZ can be expressed with the help of lines of the similar values, for instance. In the Figure 2 
we show only two maximums of the p mm corresponding to modal values of bimodal distribution of 
velocity and the maximum of zenith radiant angle value z R o = 45°. Figure 2a corresponds to mo = 10' 12 g, 
Figure 2b to mo = 10' 9 g. The regions of aerosols and meteors are separated by solid line according to 
inequality (13). 



Figure 2. The separation of meteoroids onto aerosols and meteors. Letter A and dark gray color 
corresponds to aerosols, otherwise to meteors (letter M and light gray color). The picture a 
corresponds to mass m 0 = 10" 12 g, b to mo = 10' 9 g. Dashed lines describe the equal altitudes of 
stopping. Signs “+” show positions of modes of the distribution. 


We can see in Figure 2 that the region of aerosols (denoted in the picture with the letter A and 
dark gray color) is increasing while the mass is decreasing (the region of meteors denoted as M with 
lighter gray color is decreasing accordingly). Therefore, there must be a mass value lower of which all 
particles remain aerosols. The probability for a meteoroid to become a meteor is proportional to its 
velocity and cosine of zenith radiant angle. Setting according parameters to their maximal values vo = 
72.5 km/s and cos z R o = 1 we get the critical mass value mo cr ~ 1.7 X 10' 14 g. Finally, all space particles 
entering into the Earth atmosphere remain aerosols if their masses are lower than the critical value, then 
the rate of aerosols is decreasing almost down to zero while the mass is increasing up to the value of 
approximately mo = 10' 8 g. This rate g/mo) can be easily calculated with the help of the formula 

tt/2 u q( z ro) 

g A (™ o)= \ P : ( z R 0 ) dz R 0 | P,A u o) du o • 

0 U 0MIN 


188 


The power coefficient in cumulative mass distribution k from (16) is running a range of values: k = 
0.892 for m 0 < 1.7xl0' 14 g, then k = 1.087 for m () ~ 10' 13 / 10' 12 g ,k= 1.189 for m 0 ~ 10' 11 / 10 10 g, k = 
1.438 for mo ~ 10' 9 / 10‘ 8 g (the average value for all the mass range mo ~ 10‘ 13 / 10" 8 g is k= 1.232). 

3.3 Transformation of the Primary Distribution to New Variables 

The relation (14) connects four variation parameters: three ones included into the primary distribution, 
and the fourth one, the altitude of stopping Hs. This parameter can be called conditionally the “free” 
parameter. The dotted curves corresponding to some of its values (minimal altitude of stopping, minimal 
and maximal velocities altitudes) to be expressed in kilometers are shown in Figure 2. 

Since the main aim of our investigations is to plot the two-dimensional distribution p m ii(mo, Hs) 
of the aerosol formation into the atmosphere, we will solve it in two steps. The first one is to change the 
variable z R o in the primary distribution p, m ,u(m(h uo, z ro) to the altitude of stopping H s of the aerosol 
particle. The second step will be consisting in the reducing of 3 -dimensional distribution p mv /i{mo, vo, Hs) 
to pmdrno, Hs) by means of integration of p m vn(mo, vo, Hs) over all range of v i} . 

According to statistical probability density distribution transformations and taking into account 
that only one variable is changing (zro—* Hs) we can write 


Ph(H s ) P-(H s ) ( Z R0 (Hs )) 


(H s ) 

dH s 


( 22 ) 


where zro(Hs) and determinant of transition 


d z R 0 ( H s) 

dH s 


can be found from (14). Let us denote 


C z (m 0 , u 0 , H s ) — cos z R0 C v p A (H s )/ m 0 In v 0 jv s . 


Then we can write 

p z (z R0 (H s )) = 2C z J(l-C 2 z ) , -^ = —~i=L== 
and put them into (22): 

2 2 

Ph ( m o j Li j V.s ) ~ , . * t7 z (m 0 , o 0 , H s ) 


(23). 


The final view of the obtained distribution p mv iim 0 , u 0 , Hs) while taking into account (23) is 
shown in Figure 3 for the same masses as in Figure 2. The “free” parameter now is the cosine of the 
zenith radiant angle, and the dashed curves correspond to different values of z R0 expressed in degrees. 
The value zro = 0 is placed lower than others in Figure 3 and shown with a solid curve. The region to be 
placed lower than value zro = 0 is forbidden for both aerosol and meteor particles. 

An interesting fact is that the inequality (13) is now transformed in the “stable” state in new 
coordinates and does not depend on the mass: 


H s > H* In 


CyP A m>l 

v C T bi(t>o / t>s ) y 


It is shown in Figure 3 with a solid diagonal line. 
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Figure 3. The same distributions as in Fig. 2 but in new variables. Dashed lines describe the equal 
radiant zenith angles of meteoroids. 


3.4 Resultant Distribution Reducing 

The final transformation of the distribution (23) consists in reducing it to the two-dimensional state by 
means of integration over vq. The limits of the integration can be easily determined from Figure 3 and 
according formulae. We obtain the following 


u 2 (m 0 ,H s ) 

P„A m ^ H s) = Pn, K) | Pu (°o)P H ( m 0 > °0 ’ H S )^°0 ' 


(24) 


If we denote the integral, which has to be taken numerically, as 


Dq2 ( m o s ) 




, „ , ln 2 (u n / D c ) 

^oi Oo ^ 0 SJ 


the final formula for the formation of the aerosol of space origin in the atmosphere can be written as 


P mH ( m ^ H s) = P,„ K)ttv 

ti 


C v p {) cxp(— /-/ s ./// ' ) 


m n 


1/3 


Um„H s ). (25) 


The function p m j^mo, Hs) is shown in Figure 4. 
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1E+5 p(m,H) 



Figure 4. Final two-dimensional distribution of influx of aerosol of space origin into Earth’s atmosphere. 10- 
logarithm scale 


4 Conclusion 

As it can be seen from the Figure 4 the minimal altitude which can be reachable by the aerosol stone 
particle of space origin is approximately ~ 79.6 km. This value corresponds to the meteoroid which is 
moving vertically with the velocity -16.6 km/s. A particle of the same mass and with lower velocity will 
stop higher, with higher velocity will transform into a meteor. 

The meteoroids with the mass less or equal to -1.7 x 1CT 14 g remain the aerosols always. For 
masses 10' 14 / 1CT 8 g cumulative distribution coefficient k increases from 0.892 to 1.438 while the mass 
increases. 

14 8 

The Figure 4 also demonstrates that aerosols of mass range 10' / 10' g stop in relatively thin 
altitude range 80-120 km. Evidently, the aerosols do not stay at these altitudes forever but immediately 
start to move downwards under gravitational force and the resistance force of air, which can be 
described by Stokes formula. Flow it occurs is the goal for the future work. 
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